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A  chemical  oxidation  method  was  employed  to  measure  the  kinetics  of  lithium  release  from  LiFeP04 
during  oxidation.  Similar  to  potential  step  measurements,  the  chemical  method  simplifies  quantification 
compared  to  the  common  electrochemical  techniques  (PITT,  GITT  etc.).  It  was  found  that  the  overall 
release  of  lithium  fits  one  dimensional  diffusion  kinetics,  however,  it  is  also  shown  that  the  mechanism 
must  be  more  complex  as  the  derived  activation  energy  led  to  an  unusually  low  attack  rate  of  ~  108  Hz.  A 
comparison  of  carbon  coated/carbon  free  LiFeP04  samples  indicated  that  the  carbon  coating  has  only  a 
marginal  effect  on  the  delithiation  kinetics. 

©  2014  Published  by  Elsevier  B.V. 


1.  Introduction 

One  of  the  serious  issues  facing  widespread  use  of  electric  cars  is 
the  battery  charging  time  [  1  ] .  A  large  fraction  of  lithium-ion  battery 
research  is  therefore  focused  on  understanding  kinetics  of  the 
processes  that  limit  rate  performance  of  Li-ion  batteries  [2],  This  is 
particularly  true  for  lithium  iron  phosphate  (LiFeP04),  which  is 
targeted  for  electric  vehicle  applications  due  to  its  low  price,  high 
safety,  and  above  average  practical  energy  density  ~520  Wh  kg  1 
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[3] ,  This  material  has  low  intrinsic  electronic  and  ionic  conductiv¬ 
ities  [4],  which  can  be  overcome  by  adding  conductive  coatings  to 
the  particle  surface  [5],  and  by  reducing  the  lithium  ion  transport 
path  inside  the  active  material  using  such  techniques  as  nanosizing 
the  particles  [6-8],  changing  the  LiFeP04  morphology  to  porous 
microspheres  [9],  and  by  forming  nanowires  [10],  Yet,  in  spite  of 
extensive  research,  as  well  as,  the  commercialization  of  LiFeP04,  its 
redox  mechanism  is  still  a  subject  of  intense  debate  [11],  E.g.  the 
mechanism  of  the  phase  change  between  olivine  -  LiFeP04  and 
heterosite  -  FeP04  inherent  to  the  redox  process  has  been  explained 
by  several  different  models  beginning  with  the  core— shell  model 

[4] ,  the  shrinking-core  model  [12],  and  the  more  recent  domino- 
cascade  model  [13]  etc.  Experimental  evidence  for  these  models 
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rely  on  X-ray  diffraction  (XRD)  [14,15],  Mossbauer  spectrometry 
[15,16],  X-ray  Photoelectron  Spectroscopy  (XPS)  [17],  Fourier 
Transform  Infrared  Spectroscopy  (FTIR)  [18],  Time-of-Flight  Sec¬ 
ondary  Ion  Mass  Spectroscopy  (TOF-SIMS)  [19],  Electron  Energy 
Loss  Spectroscopy  (EELS)  [20],  optical  microscopy  [20],  and  Scan¬ 
ning  Transmission  Electron  Microscopy  (STEM)  [21],  These  mostly 
ex-situ  techniques  have  been  supplemented  with  electrochemical 
measurements,  including  galvanostatic  intermittent  titration 
technique  (GITT),  electrochemical  impedance  spectroscopy  (EIS), 
potentiostatic  intermittent  titration  technique  (PITT)  and  cyclic 
voltammetry  (CV)  to  obtain  kinetic  information.  These  electrody¬ 
namic  measurements  are  intrinsically  difficult  to  perform  reliably. 
This  is  partly  due  to  the  multi-component  structure  of  the  working 
electrode,  which  consist  of  current  collector,  binder,  conductive 
carbon,  and  active  material  [22],  As  such,  the  porosity,  thickness, 
and  connectivity  can  significantly  influence  the  electrochemical 
performance  and  thus  the  measurement  [23],  In  fact,  it  is  difficult  to 
justify  the  assumption  of  unrestricted  mass  transport  in  porous 
electrodes  where  electrolyte  transport  paths  are  tortuous  [24],  The 
LiFePC>4/FeP04  system,  is  further  complicated  by  the  biphasic  na¬ 
ture  of  the  redox  reaction,  which  leads  to  a  flat  potential  plateau 
[12],  In  concert  with  oversimplified  models  [25]  a  wide  range  of 
lithium  diffusion  coefficient  (10-12  to  10"18  cm2  s  ',  see  Table  SI) 
have  been  reported  [25-29],  The  lack  of  consistency  in  literature 
makes  quantitative  prediction  of  the  battery  performance  prob¬ 
lematic  [30—35],  Further,  LiFePCU  requires  conductive  coating  to 
function  properly  in  lithium  ion  battery  electrodes,  yet,  knowledge 
about  how  the  coating  process  affects  the  lithium  transport  inside 
the  particle  is  still  elusive  [36],  Resolving  this  issue  will  require  a 
technique  that  is  applicable  to  both  coated  and  uncoated  samples 
as  presented  here. 

In  this  study  a  solution  based  chemical  oxidant  (H2O2)  [37]  is 
used  to  derive  the  critical  reaction  parameters  for  the  delithiation 
process,  via  Eq.  (1). 

V2  H202(aq.)  +  CH3COOH(aq.)  +  LiFeP04(s) 

^H20  +  CH3COOLi(aq.)  +  FeP04(s)  '  ‘ 

This  process  was  chosen  because  H202  is  highly  miscible  with 
water,  so  that  a  large  excess  of  oxidant  can  be  used,  also,  solubility 
(40.8  g  per  100  g  at  20  °C  [38])  and  diffusion  coefficients 
(~10"s  cm2  s"1  [39])  of  the  lithium  salt  in  water  are  remarkable. 
Thus,  as  opposed  to  using  nitronium  (NOj)  in  organic  solvents  [40— 
42],  here  most  concerns  regarding  lithium  transport  in  the  solvent 
are  eliminated.  Moreover,  the  extracted  lithium  can  be  directly 


quantified  using  a  standard  analytical  technique  (e.g.  atomic 
emission)  by  removing  aliquots  at  fixed  intervals  during  the 
chemical  oxidation  reaction,  which  can  further  be  correlated  to  the 
initial  mass  of  LiFePCU  [43],  This  approach  is  further  advantageous 
as  the  entire  surface  of  the  LiFePC>4  particles  are  accessed  by  the 
oxidant,  thus  minimizing  the  issue  of  non-uniform  current  distri¬ 
bution  and  shadow  effects  inherent  to  standard  electrochemical 
techniques  with  composite  electrodes. 

2.  Material  and  methods 

2  g  of  LiFePC>4  was  added  to  300  mL  of  water  containing  glacial 
acetic  acid  (10  mL,  Alfa  Aesar)  and  hydrogen  peroxide  (10  mL,  ACS 
Grade,  29.0—32.0%,  EMD).  The  suspension  was  vigorously  stirred 
and  small  aliquots  were  periodically  removed  through  a  filter  to 
immediately  separate  the  particles  from  the  reaction  medium.  This 
process  continued  until  1800  s.  The  concentration  of  lithium  in  the 
aliquots  was  analysed  by  atomic  emission  (7  =  670.8  nm)  using  a 
Varian  AA-1475  and  a  calibration  curve  (0.1 -0.8  ppm)  prepared 
from  a  standard  solution  (Specpure®,  Li  1000  pg  ml"1,  Alfa  Aesar). 
After  completion  of  the  reaction,  the  suspension  was  filtered  and 
rinsed  with  3  x  30  mL  of  water.  The  delithiated  LiFeP04  was  sub¬ 
sequently  dried  at  60  °C  under  vacuum  before  further  XRD  and 
electrochemical  analysis  (see  Supporting  information).  The  same 
procedure  was  used  for  all  LiFePC>4  samples.  The  LiFePCU  samples 
were  donated  by  Clariant  (Canada)  inc.,  Saint-Bruno  de  Montarville, 
Canada.  The  LiFePC>4  samples  were  used  as  received. 

3.  Results  and  discussion 

The  delithiation  experiments  were  performed  on  suspensions  of 
three  types  of  LiFePCU  with  different  average  particle  sizes  (see 
Supporting  information,  S3).  Type  A  (davg  =  591  nm)  is  carbon 
coated  and  formed  by  a  solid  state  reaction,  while  type  B 
(daVg  =  243  nm)  and  type  C  (daVg  =  313  nm)  are  formed  by  hydro- 
thermal  synthesis,  with  and  without  carbon  coating,  respectively. 

3.1.  Delithiation  behaviour 

LiFeP04  was  removed  from  the  oxidation  solution  by  filtration 
following  the  kinetic  delithiation  experiments,  and  subjected  to  X- 
ray  diffraction  (XRD)  (Figs.  SI  and  S2),  and  electrochemical  char¬ 
acterization  (carbon  coated  types  only).  XRD  confirmed  complete 
transformation  to  heterosite  FePCTj,  where  no  additional  peaks  were 
observed.  Iron  atomic  absorption  spectroscopy  analysis  of  the 


Fig.  1.  Deinsertion  of  lithium  C-LiFeP04  solid  state  (red  square),  hydrothermal  uncoated  LiFeP04  (black  triangle),  and  hydrothermal  C-LiFeP04  (blue  square).  (For  interpretation  of 
the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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oxidation  solution  was  used  to  confirm  that  lithium  exclusively  was 
dissolved  during  the  oxidation  process,  and  standard  coin-cell  tests 
confirmed  that  the  chemical  oxidation  process  had  no  adverse  ef¬ 
fect  on  the  materials  electrochemical  insertion/deinsertion  per¬ 
formance  (Fig.  S3).  Preliminary  analysis  of  the  kinetic  data  (Fig.  1) 
indicate  complete  delithiation  after  400  s,  350  s,  and  180  s  of 
chemical  oxidation  for  type  A,  B,  and  C,  respectively.  These  values 
did  not  change  when  the  H202  concentration  was  increased  by  a 
factor  of  10,  indicating  that  the  reaction  was  not  limited  by  trans¬ 
port  of  the  oxidant  to  the  particles  (see  Table  S2).  Further,  the 
complete  delithiation  of  LiFePC>4  (x  ~  1  in  Li(i_X)FePC>4)  indicate  a 
low  concentration  of  antisite  defects  in  the  starting  materials 
[44,45], 

3.2.  Kinetic  analysis  of  the  delithiation  reaction 


To  provide  a  more  detailed  analysis,  the  kinetic  data  were  fitted 
to  a  one  dimensional  diffusion  model  (see  Supporting  information), 
i.e.  no  nucleation  limitations  consistent  with  single  phase  deli¬ 
thiation  process  along  the  diffusionally  fast  fa-axis  [46],  This  follows 
the  Delmas  et  al.  proposed  model  where  the  chemical  delithiation 
is  initiated  at  multiple  points  throughout  the  crystallites,  i.e. 
instantaneous  nucleation  [13],  driven  by  an  overpotential  much 
larger  than  the  potential  barrier  for  the  formation  of  a  single  so¬ 
lution  region.  Malik  et  al.,  in  a  purely  computational  study,  found 
this  barrier  to  be  about  30  mV  [47],  Experimentally,  Oyama  et  al., 
based  on  potential  step  discharge  curves,  found  that  the  delithia¬ 
tion  reaction  proceeds  via  a  one  phase  process  when  a  large  over¬ 
potential  (~150  mV)  was  applied  [48],  The  overpotential  for  the 
H202  process  used  here  is  estimated  to  be  ~  1.3  V,  by  comparing  the 
potential  of  the  half-reaction  at  standard  condition. 

Using  the  single-phase  model  and  solving  Fick’s  2nd  law  in  one 
dimension  Eq.  (2)  is  obtained: 


(2) 


where  a  is  the  sample  conversion  fraction,  D  the  diffusion  coeffi¬ 
cient,  R  the  particle  radius,  and  t  the  time.  Within  this  framework 
an  apparent  diffusion  coefficient  can  therefore  be  found  from  the 
slope  of  the  ln(l  -  a)  vs.  t  plots  (Fig.  2). 


(3) 


square),  hydrothermal  uncoated  LiFeP04  (black  triangle)  and  hydrothermal  C— LiFeP04 
(blue  square).  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 


In  this  way,  apparent  one-dimensional  diffusion  coefficients  of 
2.9  x  10-12,  1.2  x  1(T12,  and  (1.4  ±  0.6)  x  10-12  cm2  s-1  (95% 
confidence  limit)  were  calculated  for  sample  types  A,  B,  and  C, 
respectively.  For  comparison,  these  values  can  be  further  converted 
to  three  dimensional  radial  diffusion  coefficients  using  approxi¬ 
mation  D3D  ~0.3  Did,  as  proposed  by  Thorat  et  al.  [30],  On  this 
basis,  diffusion  coefficients  estimated  for  the  types  B  and  C, 
exhibited  no  significant  difference.  Since  these  samples  are  essen¬ 
tially  carbon  coated  and  carbon  free  versions  of  the  same  product, 
the  presence  of  carbon  coating  only  marginally  affects  the  kinetics 
of  the  tested  samples  (Table  1 ). 

The  thermal  activation  of  the  chemical  delithiation  was  also 
examined  at  various  temperatures  (Fig.  3)  and  the  corresponding 
linear  fit  (first  order)  profiles  are  presented  in  Fig.  4. 

This  yielded  an  Arrhenius  diffusional  activation  energy  (Ea)  of  ca. 
18.4  kj  mol-1  (Fig.  5),  which  is  lower  than  the  most  cited  values  (see 
Table  SI ).  For  example,  the  activation  energy  (39  kj  mol  1 )  obtained 
by  Takahashi  et  al.  using  cyclic  voltammetry  [49]  suggests  that  the 
thermal  activation  of  phenomena  outside  the  active  material  par¬ 
ticles  played  a  role  in  their  study.  These  could  include  the  con¬ 
ductivity  of  lithium  into  the  electrolyte  [50]  as  well  as  the 
interconnection  between  particle  [51].  Our  calculated  activation 
energy  falls  within  the  range  proposed  by  Maxisch  et  al.  who  used 
Ab  initio  calculations  (17—20  kj  mol-1)  [52], 

Values  of  the  activation  energy  and  the  diffusion  coefficient  can 
be  used  to  calculate  the  product  of  the  hopping  length  (a)  and  the  v* 
attempt  frequency  (Hz),  Eq.  (4). 
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D  =  a2v*  exp(-E3/kBT)  (4) 

For  systems  with  slow  surface  reaction  great  care  should  be 
taken  when  interpreting  mechanistically  the  experimental  values 
of  the  pre-exponential  and  the  activation  energy  [53,54],  Assuming 
hopping  lengths  of  3-10  A,  which  ensures  that  all  neighbouring 
lithium  sites  can  be  reached,  an  attack  rate  of  v*  ~  108  Hz  would  be 
obtained.  This  value  is  about  5  orders  of  magnitude  smaller  than 
the  one  expected  for  a  diffusion  limited  process.  Therefore,  to  better 
represent  the  system,  a  contributory  factor  is  frequently  added  to 
the  equations,  which  represents  the  role  of  a  second  phase,  charge 
transfer,  structural  stresses,  or  a  heterogeneous  surface  energy  [53], 
We  therefore  believe  that  the  one  dimensional  diffusion  model  is 
too  simplistic  to  describe  the  mechanism  of  the  delithiation  reac¬ 
tion.  Nevertheless,  we  provide  an  apparent  diffusion  coefficient 
here  for  comparative  purposes,  showing  that  it’s  numerical  value  is 
significantly  higher  than  expected  [55]. 

3.3.  Data  analysis  using  the  standard  models  frameworks 

A  number  of  previous  kinetic  studies  have  applied  the  classical 
solid-state  kinetics  that  are  based  on  closed  form  solutions  to 
simplistic  growth  models  [56].  As  such,  Allen  et  al.  [57]  used  the 
Avrami— Johnson— Mehl— Eroofev  (JMAEK)  approximation,  Eq.  (5). 

-ln(l  -  a)1/"  =  (kt)  (5) 

where  a  is  the  sample  conversion  fraction  (here  equivalent  to  x),  k  is 
the  rate  constant,  and  n  is  the  Avrami  exponent  [53],  On  this  basis, 
Allen  et  al.  [57]  and  Oyama  et  al.  [48]  found  n  ~  1,  which,  within  the 
JMAEK  framework,  is  consistent  with  the  instantaneous  nucleation 
combined  with  a  parabolic  (diffusional)  growth  law  in  two  di¬ 
mensions  or  one  dimensional  linear  growth  [53],  The  latter  case  is 
further  consistent  with  the  moving  phase  boundary  model  as  found 
in  the  domino-cascade  model  of  Delmas  and  coworkers  [13],  Eq.  (2) 
with  n  =  1  is  also  equal  to  the  solution  to  Fields  2nd  law  in  one 
dimension  as  detailed  earlier.  Consequently,  the  kinetics  data 


Li  diffusion  in  LiFeP04  and  the  size  of  the  particle. 

Source  LiFeP04  Average  diameter  Coefficient  of  diffusion 

(nm)  10-12  (cm2  s-*) 

Solid-state  591  2.9 

Hydrothermal  (uncoated)  313  1.2 

Hydrothermal  (coated)  243  1.4 


Table  2 

Rate  equations  for  solid  phase  reactions. 


Name 

Mathematic  function  integrated  g(x)  =  kt 

Abrr. 

Avrami-Erofeyev 

[-ln(l  -  x)]»'2> 

A2 

Avrami-Erofeyev 

[— ln(l  -  x)](1,3) 

A3 

Avrami-Erofeyev 

[— ln(l  -  x)](I,4) 

A4 

ID  diffusion 

X2 

D1 

2D  diffusion 

((1  —  x)ln(l  -x))  +  x 

D2 

3D  diffusion  (Jander) 

(l-0-x)W3>)2 

D3 

Ginstling-Brounshtein 

1  -  (2/3  )x  -  (1  -  x)(2,3> 

D4 

Contracting  area 

1  -(1  -x)«2> 

R2 

Contracting  volume 

l-d-x)'1® 

R3 

P2 

P3 

Power  law 

x1/4 

P4 

First-order 

-ln(l  -  x) 

FI 

Second-order 

[1/(1  -  x)]  -  1 

F2 

Third-order 

(1/2)[(1  -x)-2  -  1] 

F3 

cannot  be  used  to  confirm  a  moving  barrier  or  a  single  phase  re¬ 
action  mechanism  without  additional  information.  For  complete¬ 
ness,  we  also  tested  the  standard  steady  state  diffusional  growth 
models  (Table  2;  D1-D4). 

Importantly,  we  do  not  find  a  good  fit  for  a  ID  moving  reaction 
front  coupled  with  diffusion  model,  which  would  be  a  consequence 
of  Li+  movement  predominately  along  the  b-axis  channels  after 
having  been  released  at  the  LiFePCU/FePCU  phase  boundary.  How¬ 
ever,  fitting  the  time  resolved  data  to  a  3D  moving  reaction  front 
with  diffusion  model  yielded  fits  of  similar  quality  as  those  found 
for  Eq.  (2)  (Table  3;  D3,  D4).  This  type  of  model  has  been  success¬ 
fully  invoked  in  describing  the  early  stages  of  constant  current 
discharge  curves,  which  was  rationalized  in  terms  of  a  reaction 
front  at  the  agglomerate  scale  by  Zhu  et  al.  [25],  Further,  2D  moving 
reaction  front  with  diffusion  model  is  indistinguishable  from  the 
3D  model  within  the  presented  experimental  data. 

4.  Conclusion 

The  chemical  oxidation  for  kinetic  studies  of  delithiation  pre¬ 
sented  in  this  report  benefits  from  an  instantaneous  potential  step 
to  the  entire  sample  while  allowing  unrestricted  transport  of 
lithium  ions  from  the  oxidized  particles.  As  such,  it  overcomes  two 
of  the  most  important  issues  associated  with  the  standard  GITT/ 
PITT  technique,  i.e.  a)  the  potential/current  is  not  uniform  due  to  the 
electrolyte/electrode  resistance  and  b)  ionic  transport  is  hampered 
due  to  the  tortuous  path  of  the  electrolyte  in  the  composite  elec¬ 
trode.  Consequently,  the  apparent  diffusion  coefficients  found  in 
this  study  are  comparable  to  the  highest  values  reported  in  litera¬ 
ture.  That  is,  our  results  were  even  higher  than  the  values  reported 
(1.9  x  10-13  cm2  s_1)  by  Yu  et  al.  who  used  an  NOj  in  acetonitrile 
oxidant  flow  through  a  packed  bed  of  LiFeP04  to  drive  the  deli¬ 
thiation  reaction  [40],  This  suggests,  surprisingly,  that  the  liquid 
flow  was  insufficient  to  support  unobstructed  oxidant  and/or 
lithium  transport  between  the  liquid  and  solid  phase.  Further,  we 
provide  here  for  the  first  time  uncoated  LiFePCU  oxidation  data, 
where  complete  delithiation  takes  place  within  ~400  s  for 
~300  nm  diameter  particles.  This  indicates  that  the  role  of  state- 
of-the-art  carbon  coating  as  an  inhibiter  of  lithium-ion  transport 
should  be  negligible.  Finally,  while  there  is  substantial  indirect  and 
theoretical  evidence  for  a  single  phase  delithiation  mechanism  at 
high  overpotentials,  a  conclusive  dynamic  structural  analysis  is  still 
missing.  This  is  likely  due  to  the  insufficient  ionic  and/or  electronic 
charge  transport  to  and  from  the  particles  during  the  reaction  using 
standard  techniques.  The  methodology  elaborated  here,  together 
with  an  appropriate  structure  analysis  technique  could  provide  a 


D.  Lepage  et  al.  /  Journal  of  Power  Sources  256  (2014)  61—65 


Table  3 

Correlation  coefficients  (r2)  for  each  rate  equation  solid  phase  reaction  (best  fit  are  highlighted  in  bold)  AVG  =  average. 


Models  Hydrothermal  coated  Hydrothermal  coated  (32  °C)  Hydrothermal  coated  (5  °C)  Solid  state  Hydrothermal  Unc.  AVG 


A2  0.937  0.912 

A3  0.956  0.937 

A4  0.987  0.948 

D1  0.977  0.980 

D2  0.995  0.990 

D3  0.999  0.987 

D4  0.998  0.991 

R2  0.974  0.972 

R3  0.984  0.982 

P2  0.902  0.888 

P3  0.891  0.869 

P4  0.884  0.858 

FI  0.996  0.994 

F2  0.976  0.971 

F3  0.945  0.875 


0.948 

0.968 

0.976 

0.977 

0.996 

0.998 

0.999 


0.891 

0.885 

0.999 

0.960 

0.868 


0.887  0.927  0.922 

0.923  0.951  0.947 

0.964  0.962  0.967 

0.961  0.984  0.976 

0.987  0.997  0.993 

0.993  0.988  0.993 

0.993  0.997  0.996 

0.957  0.976  0.972 

0.973  0.986  0.983 

0.830  0.887  0.882 

0.806  0.871  0.865 

0.794  0.862  0.856 

0.993  0.995  0.995 

0.945  0.939  0.958 

0.800  0.808  0.859 


powerful  tool  toward  understanding  the  intricacies  of  the  complex 
LiFeP04  lithium  deinsertion  process. 
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